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As one of the fitting methods, the polynomial approximation is effective to process sophisticated 
problem. In this paper, we employ this approach to handle the scattering of scalar field around 
the Schwarzschild-de Sitter black-hole. The complex relationship between tortoise coordinate and 

© : radial coordinate is replaced by the approximate polynomiaL The SchrMinger - like equation ' the 

real boundary conditions and the polynomial approximation construct a full Sturm-Liouville type 

(N 
On 



> 
(N 

m 



O 



problem. Then this boundary value problem can be solved numerically according to two limiting 
cases: the first one is the Nariai black-hole whose horizons are close to each other, the second one 
is when the horizons are widely separated. Compared with previous results (Brevik and Tian), the 
field near the event horizon and cosmological horizon can have a better description. 
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I. INTRODUCTION 



In the black-hole physics Hawking radiation [1( is always a very important conception which indicates that black- 



holes are not perfect black, but radiate thermally and eventually explode. An in-depth discussion of derivation can be 
obtained in Ref . [2] As for the recent research, one can refer to Ref. Many researchers have developed various methods 



O . 

and techniques to study the black-hole by using the radiating particles, such as the simple Klein-Gordon particles 
and Dirac particles (for some early works, see Damour and RuffiniQ] and Chandrasekhar[5| respectively). Recently, 
as for scalar field, Higuchi et al. 6] and Grispino et al.Q gave its solution outside a Schwarzschild black- hole, Brady 



any works, see Damour ai 
■ ao iui o^aiaj. uciu, in 5 uuii cu ai. [fj and Grispino et al.Q 
$_i ' et al. Q] studied the Schwarzschild-de Sitter case and Guo et al.[9( made further studies in the Reissner-Nordstrom-de 
Sitter one. 

The Schwarzschild-de Sitter (SdS) space is a spherically symmetric system. 10] It can be treated as a small 
Schwarzschild black-hole embedded in de Sitter universe. In this space there are two horizons: one is inner black-hole 
horizon r e and the other is outer cosmological horizons r c . In 2001, Brevik and Simonsen[ll] gave a massless scalar 
field solution by tangent approximation which contains an explicit tangent function. Viewing from the global frame, 
this method matches r with f very well, where r is the radial coordinate and f is a fitting function. However, its 
insufficiency is the weak fitting near the two horizons. Even in the intermediate zone, the fitting r with f is not precise 
enough. In their paper, [ill they have studied for two extreme cases whose horizons are either very close to each other 
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or lie very far away. Afterwards, Tian et al 



12| used a polynomial approximation containing 20 monomials and gave 



another different numerical solution only in the extreme Nariai black-hole. This useful polynomial approximation 
is more precise than the tangent approximation. Especially, in the leading intermediate zone, the fitting r with a 
polynomial makes a good match. However, this type of approximation rapidly deteriorates near the horizons, i.e. 
inappropriate boundary conditions are used. On the other hand, the widely separated horizons case is not considered 
in Ref . [121 ] Considering the above situation, we re-study the scattering of a scalar field with polynomial approximation. 

This paper is organized as follows: in Section 2, we present the Schwarzschild-de Sitter space and point out the 
positions of black-hole horizon and cosmological horizon. In Section 3, by the polynomial approximation, the full 
Sturm-Liouville type problems are solved for two extreme cases. Section 4 is a conclusion. We adopt the signature 
(+, — , — , — ) and put h, c ,and G equal to unity. The same setting of field parameters is employed in Ref. 11 1 



II. SCALAR FIELD IN SCHWARZSCHILD-DE SITTER SPACE 



The spherically symmetric metric of the Schwarzschild-de Sitter space [lfj is given by 

1 

lir 



ds 2 = f(r)dt 2 - -^-dr 2 - r 2 (d6 2 + sin 2 9d<f> 2 ) , (1) 



where 

/(r ,_,_2«_^, (2) 

r 3 

with the black-hole mass M and the cosmological constant A. It is an exact exterior solution of the Einstein field 
equations for a spherical mass distribution 

R^v - ^SV-R + As^ = (3) 

Here, we take the cosmological constant A as a free parameter. The similar process can be found in Refs. [ill H. Il^ 
This space is bounded by two horizons — an inner horizon (black-hole horizon) and an outer horizon (cosmological 
horizon). Under the limit A — ► 0, this metric has exactly the same line-element as the Schwarzschild space. But in 
the limit of M — > 0, it reduces to the de Sitter one. 
Mathematically, expression @ can be rewritten as 

f(r) = —(r-r e )(r c -r)(r-r ). (4) 

The singularity of metric JT]) is determined by f(r) = 0. The solutions to this equation are these to inner horizon r e 
and outer horizon r c , as well as a negative solution r a = — (r e + r c ). The last one has no physical meaning. Here we 
only consider the positive solutions. The positions of r c and r e are given by 

rc = _^ cos??; 

e = ^COs(120°-7 ? ), 

where r\ = 1/3 arccos(— 3My/X) with 30° < r\ < 60°. The real physical solutions are accepted only if A satisfies 
AM 2 < 1/9. 14( If the cosmological constant A reaches its maximum, the Nariai black-hole appears. 
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A massless scalar field r, 0, 0) is considered here. Using the separable solutions[l 

$ 

the scalar field equation 

□$ = 0, 

is decomposed into two differential equations: 



(6) 
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where Ruj(r,t) is the time-dependent radial function and Yi u (9,(f>) is the spherical harmonics function. Equation 
(jSJ) determines the evolution of evaporating black-hole. It is necessary to eliminate the time variable by the Fourier 



component e via 



So Eq.© can be rewritten as 



Rw(r,t) -> *^(r)e" 



-f{r)Uf(r)^) + V(r) 
ar dr 



Vu,l(r) = w 2 ^/(r), 



whose potential function is given by 

V(r) = f(r) 
Now we introduce the tortoise coordinate 



ldf(r) | 1(1 + 1) 
r dr r 2 



1 



dr 



2MJ f(rY 

The tortoise coordinate can be expressed by surface gravity as follows: 



(10) 
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Explicitly, we have 
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FIG. 1: The radial coordinate r (full line) and polynomial approximation y (triangle-up) versus the tortoise coordinate x with 
A = 0.11. 



By tortoise coordinate transformation (jT2J) , the radial equation (fTTj) can be written in the "Regge - Wheeler" form 

* B ,W=4MV* ul (4 (19) 



which has the form of Schrddinger equation of quantum mechanics. It is usually called Schrodinger-like equation. 
The incoming or outgoing particle flow between inner horizon r e and outer horizon r c is reflected and transmitted by 
the potential barrier V(r). The evolution of wave solution of massless scalar field is also determined by potential 
V(r). 

III. FULL STURM-LIOUVILLE TYPE PROBLEM 

According to Eqs. (|1]) and (jT2J) , the potential V(r) disappears near black-hole horizon r e and cosmological horizon 
r c . So the potentials near horizons are given as follows: 



V(r e ) = V(r c ) = 0. 



Hence near the horizons, Eq. (JT^J) reduces to 



dx 2 



4M 2 w 2 



* W i(a?) = 0. 



Absolutely, its solutions are e ± l2MuJX or their comprehensive form. Taking the real scalar field[ll|, [l^l 
we choose the real part of its solutions as the boundary condition 



cqs{2Mujx) 



(20) 
(21) 

into account, 
(22) 



There are two coordinates — radial coordinate r and tortoise coordinate x contained in Eq. (|19[) . However, the 
source transformation expression (114|) is too complicated to invert it to the form of r — r(x). In order to solve Eq. (fl9|) 
conveniently, it is necessary to use an approximate method for transition 



f(x), 



(23) 
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FIG. 2: The potential V(x) versus tortoise coordinate x with A = 0.11, M = 1, I = 1, where the maximum value is represented 
by the black point. 



where the radial coordinate r is replaced by the fitting function f containing the tortoise coordinate x. Here we use 
the polynomial approximation [2] to fit r with f . This method involves complicated polynomials unlike the explicit 
tangent approximation. [llj For any given value of the parameter A, we can always find an appropriate approximate 
method from the above by adjusting the parameters. The difficulty is how to obtain the fitting function (|23|) . Here, 
we consider the polynomial approximation 12 



JV 

y = f = ^2 QiX% 

i=0 



(24) 



where is the coefficient and N is the degree of polynomial. In the various fitting functions, the greatest advantage 
of the polynomial approximation is that it can obtain the optimal approximation by the adjustable parameter N. 
One should note that it is wrong to use the bigger N to gain a more accuracy approximation. N must be chosen 
according to the fitting interval. By combining the potential V(x), the Schrodingcr-like equation (fl"9]l . fitting function 
(|24|l with the boundary conditions (126j) , we can present a full Sturm-Liouville type problem. This kind of boundary 



value problem is usually used to solve the field equation, such as Refs.flll [l2. Il3l| Note that the same setting of field 
parameters as that in Rcf.[ll| are adopted in the subsections. 



A. More Exact Boundary Conditions for Nariai Case: A = 0.11 



Nariai solution has been discovered by Nariai, which is the exact solution to the Einstein equation with A > 
without a Maxwell field. Its topological structure is a (l+l)-dimensional dS spacetime with a round 2-sphere of fixed 



11 



radius, i.e., (IS2 x S 2 . We adopt the same value A = 0.11 as appeared in Refs 
we find the inner horizon is r e — 2.8391A/ and the outer horizon is r c = 3.1878M. Just as in Ref. 
same polynomial approximation here. The coefficients {di} are showed in Table 1. 



12j Then substituting A into Eq.flD, 



12| , we employ the 
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FIG. 3: The wave function ip^i versus the tortoise coordinate x with A = 10 , M = 1 and 1 = 1. 
TABLE I: The coefficients of degree 20 polynomial with A = 0.11. 



a = 2.9817 


oi = 6.5107 x 10" 3 


a 2 = 4.0912 x 10 -5 


o 3 = -2.9913 x 10" 6 


a 4 = -3.4895 x 10" 8 


a 5 = 1.6009 x 10" 9 


a 6 = 2.3413 x 10 -11 


a 7 = -8.0083 x 10~ 13 


a s = -1.1964 x 10~ 14 


a 9 = 3.3845 x 10~ 16 


oio = 4.3110 x 10~ 18 


a n = -1.0899 x 10" 19 


aw = -1.0120 x 10~ 21 


013 = 2.4364 x 10" 23 


oi4 = 1.4031 x 10~ 25 


oi6 = -3.4329 x 10" 27 


oie = -9.5242 x 10~ 30 


Oi7 = 2.6439 x 10 -31 


ois = 1.5652 x 10" 34 


019 = -8.0402 x 10" 36 


a 20 = 4.3262 x 10" 39 



The boundary conditions in Ref. 12] came directly from the original work in Ref. 11| . 

# w /0)U=-ioo = *u,/(aOU=ioo = cos(200Mo;). 



(25) 



However, it is not appropriate to use Eq. (|25p directly as the boundary conditions do not consider the new approxi- 
mation. The intervals of the boundary conditions should be in accord with the fitting intervals. Now, we present the 
inappropriate boundary conditions in previous polynomial approximation, [l^ One can treat the polynomial (|24p as 
a function varying with x for N = 20. The coefficients are listed in Table. 1. The functional images of Eqs. (|M|) and 
(|14|) both are drawn in Fig.l. It is shown that f (or y) and r match exactly in the intermediate zone. However, near 
the horizons r e and r c , the polynomial behavior takes over and the approximation quickly deteriorates. Especially, 
there are significant differences between r and y in the two intervals, [—70, —100] and [70, 100], along the horizontal 
axis. Because of the unnecessary intervals, the singular peak has arisen in the waves (one can refer to Figs. 4 and 5 in 
Ref. 



121]). This problem can be solved by reducing the interval x from [—100, 100] to [—70, 70]. Hence, the effective 
interval in the radial direction has been changed from [2.8391M, 3, 1878 A/] to [2.8382M, 3.1865M]. After removing 
the useless intervals: [—70, —100] and [70, 100], we obtain another exact boundary conditions 



*^(z)U=-7o = *^(z)U=70 = cos(140Alw). 



(26) 



The potential V(x) of Nariai black-hole is plotted in the range of —70 < x < 70 in Fig. 2, with maximum value 
V = 7.4381 x 10~ 4 corresponding to x = 2.0352. By using Mathematica software in book, 17j one can solve it 



numerically as a boundary value problem, where the command NDSolve is used. The amplitude versus the tortoise 
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FIG. 5: The radial coordinate r (full line), the 10-th polynomial i/io (triangle- up) and the 20-th polynomial j/20 (star) versus 
the tortoise coordinate x with A = 10 -3 . 



coordinate is shown in Fig. 3. It is seen that the solution ^^(x) is similar to a harmonic wave without considering 
the decay factor 1/r in the ansatz ©. The radial equation ([5]) is transformed into a standard wave equation by 
useful tortoise transformation (fT3|) . With the real boundary conditions (cosine functions), the harmonic wave arises 
naturally. Taking into account the actual case, we also plot the amplitude versus r in Fig. 4. This diagram illustrates 
clearly that waves stack up near r e and r c . 

For the case of Brevik's work, because the tangent approximation does not work well near the two horizons (see Fig. 3 
in Ref. 11]), waves do not pile up near the outer horizon r c (see Fig. 6 in Ref.[ll |). For the case of Tian's work, there 



is a singular peak in waves near x — 2.8582 (or r ~ 3) (see Figs. 4 and 5 in Ref. [121). The occurrence of this singular 
peak is due to inappropriate boundary conditions (|25|) chosen. In this paper, the precise polynomial approximation is 
kept, but the boundary conditions (|25p are replaced by the new exact ones (I26[) . However, viewing from the numerical 
solutions shown in Figs. 3 and 4, we can say that the afore-mentioned deficiencies have been remedied. 
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FIG. 6: The potential V(x) versus tortoise coordinate x with A = 10 3 , M = 1, I = 1. 
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FIG. 7: The variation of the wave function with the tortoise coordinate x for A = 10 , M = 1 and Z = 1. 



B. Widely Separated Horizons Case: A = 0.001 



The accelerating universe phenomenon can be easily explained by a repulsive force produced by non-zero and positive 
cosmological constant Ao ~ 10~ 52 to _2 .|18[ The cosmological constant has very interesting grav itational effects on 
various astrophysical scales such as the gravitational lensing statistics of extragalactic surveys, I19l large-scale velocity 



22j and Solar 23]). Considering the 



flows, 20] the effects on observation in small system (Galactic, j^lj Planetary, 
widely separated horizons case, we take the cosmological constant A = 1CP 3 , which is employed as a general setting 
in many works. 

QQ of course, the other value of cosmological constant subjected to the condition AM 2 < 1/9 also 



11] 



can be adopted in principle. Here we take the same value as that in Rcf. 

In the polynomial approximation, the number of terms can be selected at random in principle. Two points determine 
one line or an approximating polynomial of degree 1; three points determine an approximating polynomial of degree 
2, and so on; while n + 1 points determine the approximating polynomial of degree n. However, the polynomial with 
higher degree presents some defective numerical characteristics. With increasing degree, the fitting curve becomes 
lack of smoothness because the higher degree polynomial can be differentiated many times before it reduces to zero. 
We thus choose a right polynomial rather than a higher degree one. 

The polynomials with 20-th degree (j/20) and 10-th degree (yio) are plotted in Fig. 5. We can read that the fitting 
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TABLE II: The coefficients in the approximating polynomial of degree 10 



a = 2.3895 


a a = 0.64096 


a 2 = 0.11802 


a 3 = -6.7621 x 1(T 3 


a 4 = 1.8818 x 10~ 4 


a 5 = -3.2408 x 10" 6 


a 6 = 3.6365 x 10" 8 


a 7 = -2.6487 x 10~ 10 


a 8 = 1.2019 x 10" 12 


o 9 = -3.0754 x 10~ 15 


oio = 3.3792 x 10" 18 
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FIG. 8: The variation of the wave function ipui with the radial coordinate r for A = 10 -3 , M = 1 and / = 1. 

interval of yio is [—5, 140] and the fitting interval of j/20 is [ — 5,40]. Obviously, the former is much wider than the 
latter. Here, we adopt the 10-th degree polynomial as fitting function (f2"4"|) . The coefficients are listed in Table 2. 
Using this approximation, the potential V(x) is plotted in Fig. 6. The curve peak approaches the event horizon r e 
and is far apart from the cosmological horizon r c . Like the former case, this kind of boundary value problem can be 
solved numerically 17] with a more general setting which is more faithfully representing our world. The numerical 
solutions ^u,i(x) and ^ u i(r) are shown in Figs. 7 and 8, respectively. Obviously, with decreasing A the waves become 
much sparser near r e and much denser near r c . Otherwise, since all parts of the potential arc in the region of x > 0, 
the wave solutions ty(x) and 'F(r) are concentrated in the positive horizontal axis too. 



IV. CONCLUSION 

We have solved the real scalar field numerically with the polynomial approximation. Unlike the previous original 
work, [2 we have not only surveyed a more general setting of field parameters, i.e. the widely separated horizons 
case, but also selected more precise boundary conditions. We summarize what have been achieved as follows. 

1 Polynomial approximation is an important and comprehensible approximate method. In this paper this method 
has been used to fit the radial coordinate with the tortoise coordinate. Unlike the tangent approximation. [ill] one 
main merit of this method is that there is an adjustable parameter, the degree of the polynomial. The degree of the 
polynomial must be selected to ensure that we can obtain the maximum fitting interval. From the analysis in this 
paper, we find that the degree of the polynomial should be reduced with increasing cosmological constant. When we 
consider two extreme cases with A = 0.11 and A = 0.001, the degree should be reduced from 20 to 10. 

2 As for the Nariai black-hole, because the potential barrier V(r) (jl2[) vanishes near the two horizons, Eq. (|2"2"]) 
becomes a universal boundary condition for real solution case. But the different fitting intervals need different 
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boundary conditions. So it has to choose a well and suitable one according to the quality of fitting. In this paper, we 
have chosen another exact boundary condition. Although the previous polynomial approximation is more accurate 
to match r with f in the intermediate zone, the interval of x does not keep the original one ([—100,100]). It is 
illustrated clearly in Fig.l. Obviously, when the tortoise coordinate x is in the range of [—100, —70] or [70, 100], the 
approximation (124p quickly deteriorates. So it is necessary to use the new interval [—70, 70] to replace the previous 
one [—100, 100]. As mentioned in the above sections, by tortoise coordinate transformation the radial equation can 
be rewritten as a standard wave equation form. By combining with the real boundary condition (cosine function 
form), the harmonic wave appears. Otherwise, for the compactness of the tortoise coordinate the waves pile up near 
the horizons naturally. These effects in the Nariai case in the Brevik's work ll| are subtle and subclinical near the 
cosmological horizon r c . After rebuilding the scalar field in this paper, we find that there is no singular peak in waves, 
and the waves pile up near the two horizons r e and r c , which refreshes previous works in Refs. Ill |12 1 
3 For the widely separated horizons case, the dimensional version of A = 10~ 3 reads 

A (^) = 10-3, (2?) 

where G is the gravitational constant, c is the speed of light and M is the mass of black-hole. Considering a usual 
Stellar Black Hole, we assume the mass M equal to ten sun masses (i.e. M = 1OM ). It is known that for the Solar 
system we have 

= 1.475fcm. (28) 
Substituting the mass M and the notation (|2"8")) into Eq. (|2"Tj) , we can obtain a dimensional cosmological constant 

A = 4.6 x 1Q- I6 cm- 2 . (29) 
Using the same method, we can also obtain the dimensional cosmological constant in the Nariai case. 

A = 2.4 x IQ- Vi cm- 2 (30) 

Although the values are much larger than the observed value Ao, it is very necessary to further research them. After 
all, the space could be a Schwarzschild one if we take the observational value 10 -52 m~ 2 . 



V. ACKNOWLEDGE 



Project supported by the National Basic Research Program of China (Grant No. 2003CB716300) and National 
Natural Science Foundation of China (Grant No. 10573003). 



[1] Hawking S W 1974 Nature 248 30 

Hawking S W 1975 Commun. Math. Phys. 43 199 
[2] Brout R, Massar S, Parentani R and Spindel P 1995 Phys. Kept. 260 329 
[3] WuSQ and Cai X 2002 Chin. Phys. 11 661 

Jiang Q Q, Yang S Z and Wu S Q 2006 Chin. Phys. 15 2523 



11 

Ren J, Cao J L and Zhao Z 2006 Chin. Phys. 15 2256 

Hu Y P, Zhang J Y and Zhao Z 2007 Acta Phys. Sin. 56 683 (in Chinese) 

Jiang Q Q and Wu S Q 2006 Acta Phys. Sin. 55 4428 (in Chinese) 

Cao J L 2006 Acta Phys. Sin. 55 2682 (in Chinese) 
[4] Damour T and Ruffini R 1976 Phys. Rev. D 14 332 
[5] Chandrasekhar S 1976 Proc. Royl. Soc. Lond.A 349 571 

[6] Higuchi A, Matsas G E A and Sudarsky D 1998 Phys. Rev. D 58 104021 gr-qc/9806093 

[7] Crispino L C B, Higuchi A and Matsas G E A 2000 Class. Quant. Grav. 17 19 gr-qc/9901006 

[8] Brady P R, Chambers C M, Laarakkers W G and Poission E 1999 Phys. Rev. D 60 064003 gr-qc/9902010 

[9] Guo G H, Gui Y X and Tian J X 2003 Int. J. Mod. Phys. A. 18 4829 

[10] Rindler W 2001 Relativity (Oxford: Oxford University Press) 

[11] Brevik I and Simonsen B 2001 Gen. Rel. Grav. 33 1839 

[12] Tian J X, Gui Y X and Guo G H 2003 Gen. Rel. Grav. 35 1473 gr-qc/0304009 

[13] Liu M L, Liu H Y, Xu L X and Wesson P S 2006 Mod. Phys. Lett. A 21 2937 gr-qc/0611137 

Kanti P, Grain J and Barrau A 2005 Phys. Rev. D 71 104002 hep-th/0501148 

Liu M L, Liu H Y, Luo F and Xu L X 2007 Gen. Rel. Grav. 39 1389 gr-qc/0705.2465 
[14] Liu H Y 1991 Gen. Rel. Grav. 23 759 
[15] H Nariai 1950 Sci. Rep. Tohoku Univ. 34 160 

H Nariai 1951 Sci. Rep. Tohoku Univ. 35 62 

Liu M L, Liu H Y, Wang C X and Ping Y L 2007 Int. J Mod. Phys. A 22 4451 gr-qc/0707.0520 
[16] Jensen B P, Candelas P 1986 Phys. Rev. D 33 1590 

[17] Wolfram S 1996 The Mathematica Book 3rd ed (Cambridge: Wolfram Media/Cambridge Univ. Press) 
[18] Schmidt B P et al. [Hi-Z Supernova Team Collaboration] 1998 Astrophys. Journ. 507 46 astro-ph/9805200 

Riess A G et al.1998 Astron. J. 116 1009 astro-ph/9805200 

Peebles P J E and Ratra B 2003 Rev. Mod. Phys. 75 559 astr-ph/0207347v2 

Ostriker J P and Steinhardt P T 1995 Nature 377 600 
[19] Quast R and Helbig P 1999 Astron. Astrophys. 344 721 astro-ph/9904174 
[20] Zehavi I and Dekel A 1999 Nature 401 252 astro-ph/9904221 
[21] Whitehouse S B and Kraniotis G V 1999 astro-ph/9911485 
[22] Cardona J and Tejeiro J 1998 Ap. J. 493 52 

[23] Kagramanova V, Kunz J and LSmmerzahl C 2006 Phys. Lett. B 634 465 gr-qc/0602002v2 



